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The complete one-loop electroweak corrections to charged Higgs boson pair production in 
in ' 77 collision mode at linear colliders in the minimal supersymmetric standard model (MSSM), 

are calculated in this paper. We discuss the dependence of the corrections to the subprocess 
f~i ' 77 H^H^ on the related parameters, such as the colliding energy, charged Higgs boson 

C^, mass Mu± and some supersymmetric parameters tan /3, Msusy and gaugino mass parameter 

^ • M2. We find that the corrections generally reduce the Born cross sections and the relative 

D I one-loop corrections to both the subprocess typically in the range of —10% to —30%. We 

' also present the numerical results at the SPSla' point from the SPA project. We conclude 

that the full one-loop electroweak corrections to subprocess 77 H^H^ and the parent 
, J ■ process e^e~ — > 77 — > H'^ H~ are significant and therefore should be considered in precise 

r> \ analysis of charged Higgs boson pair productions via 77 collision at future linear colliders. 
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I Introduction 

New physics beyond the standard model (SM) has been intensively studied over the past years^. 
Most extensions of the SM require an electroweak symmetry breaking sector, which is composed 
of two scalar isospin-doublets, and the charged Higgs bosons are part of its physical spectrum 
at the weak scale. The minimal super symmetric standard model (MSSM) is one of the typical 
example. From the phenomenological point of view, if a light neutral Higgs boson was found, 
it is still very hard to tell which model it belongs to, since the fundamental properties of such 
a particle (quantum number, couplings, branching ratios, etc.) are almost the same in some 
models, e.g. in the SM and in the 'decoupling regime' of the MSSM (i.e., when Mj^o > 200 GeV, 
Mjifi ~ Mho ~ Mjj± » M/jo). While the discovery of the charged Higgs boson is an unam- 
biguous signature of existing new physics beyond the SM. 

Historically, a lot of effort has been invested in the charged Higgs boson pair production 
at the future colliders, such as the CERN Large Hadron Collider (LHC), Tevatron, and the 
proposed linear colliders (LC): NLC_2 , JLC_3 , TESLAgl and CLICg. Refs.|ni[7I presented 
the calculations of the charged Higgs boson pair productions at hadron colliders in different 
important production channels. It shows that the production cross section can reach few femto- 
bar. Linear colliders can also produce the charged Higgs pair with larger production rate, 
because the process can occur at the tree level and is not suppressed by the light Yukawa 
couplings. Furthermore, the signature of event at LC is much cleaner than that produced 
at hadron colliders. With the help of high integrated luminosity, the precise measurement 
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at LC for probing the new physics is possible. Therefore, the theoretical calculations beyond 
the tree-level are necessary in studying the charged Higgs boson productions. In Ref.[S], the 
process e~^e~ — > H^H^ involving one-loop fermion and sfermion corrections has been studied, 
it points out that the corrections are about —10% in a wide range of parameter space of the 
MSSM. Ref.[S] gives the complete one-loop electroweak corrections to the cross section of the 
process e'^e~ — > H^H^ in the THDM as well as the MSSM. It shows that the corrections vary 
in the range between —15% and 10%. The 0{a'm1 /m^f) Yukawa corrections to the process 
e^e~ — > 77 ^ H^H^ in the THDM were studied in Ref.|lUj. Ref. jllj presents the squarks 
one-loop corrections to the process e'^e~ ^ 77 — > H^H^ in the MSSM. It says that the relative 
corrections are from —25% to 25%. From the previous works which deal with the complete 
one-loop corrections to the new particle production processes, we know that the detailed study 
of the one-loop electroweak corrections for those processes at a very high colliding energy is 
necessary. An electron-positron LC can be designed to operate in either e~^e~ or 77 collision 
mode. 77 collision is achieved by using Compton backscattered photons in the scattering of 
intense laser photons on the initial polarized e~^e~ b earns Normally, the cross section for 
77 — > H^H^ is larger than that of e^e^ H^H^ due to the fact that the production rate in 
e'^e'^ collision mode is s-channel suppressed. In this paper, we present the calculations of the 
full one-loop radiative corrections to the process e'^e~ ^ 77 ^ H^H^ in the MSSM. The paper 
is organized as follows. In Sec. II. we discuss the LO results of the subprocess 77 — > H^H^ . 
In Sec. III. we give the analytical calculations of the full one-loop corrections. The numerical 
results and discussions are presented in Sec. IV. Finally, we give a short summary. 
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II The Leading Order Cross Section of subprocess 77 — ^ H^H' 



We denote the subprocess 77 — H^H as 



(2.1) 






where pi, p2 and ps, ^4 represent the four-momenta of the incoming partons and the outgoing 
particles, respectively. 

H+ 

H+ 



H- 
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Figure 1: The leading order diagrams for the 77 — H^H~ subprocess. 

The Feynman diagrams for the subprocess 77 H^H~ at the leading order(LO) are shown 
in Fig.l. There are three Feynman diagrams for this subprocess at the tree-level. The corre- 
sponding tree-level amplitude of the subprocess 77 ^ H^H~ can be represented as 



Mo = + Mo" + 



(2.2) 



where Mq, Mq and Mq represent the amplitudes arising from the t-channel, u-channel and 
quartic coupling diagrams, respectively. The explicit expressions can be written as 



(2.3) 



Mi = A4{t^u, P2^Pi), Mo« = 2ieV%(Pi)e.(P2). (2.4) 

The Mandelstam variables t, u and s are defined as t = {pi — ps)^, u = {pi — ^4)^, s = [pi +^2)^- 
Then the LO cross section for the subprocess 77 H^H^ is obtained by using the fohowing 
formula: 

a°(5,77 - H^H-) = I dt Y^IM^, (2.5) 

where imax,min = i^H± ~ 1^) ^ ^ y^s^ — 4M^±s. The summation is taken over the spins of 
initial and final states, and the bar over the summation denotes averaging over the spins of 
initial partons. 

Ill The Calculation of the Full One-loop Corrections to the sub- 
process 77 H^H 

In our calculations we use the t'Hooft-Feynman gauge. In the calculation of one-loop diagrams 
we adopt the definitions of one-loop integral functions in Ref. jl3j. In order to control the ultra- 
violet(UV) divergences, we take the dimensional reduction (DR) regularization scheme, which is 
commonly used in the calculation of the electroweak correction in the framework of the MSSM 
as it preserves supersymmetry at least at one-loop order In doing renormalization we use 
on-mass-shen(OMS) scheme [T^. The Feynman diagrams and their amplitudes are automatically 
generated by using Fey n Arts 3 package ITHl. 
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III.l Virtual Electroweak One- loop Corrections 

There are total 570 one-loop Feynman diagrams for the subprocess 77 H^H~ in the MSSM, 
and we can classify them into four groups: self-energy, vertex, box diagrams and counter-term 
diagrams. Let's consider the counter terms at first. The Higgs potential in the MSSM can be 
divided into four parts 

VH = vSUvSHvSHvi^^ (3.1) 

which represent the linear, quadratic, cube and quartic terms respectively. The linear and 
quadratic terms can be expressed as 

..(2)charge _ 1 / ^+ „-+ n t ^GG ^GA \ f '^^ \ 

where 

boG = TT^i^H cos (a - P) -Th sm{a - /?)], 

bAA = ;-[Tff (sin^ /3cosa -I- cos^ /3sina) -I- r/j(cos^ /3cosq; — sin^ /^sino;)], 

mw sin 2p 

bo A = Tr—[THSm{a-P)+ThCOs{a-P)]. (3.3) 
We use the following definitions of the renormalization constants related in our calculation 

as, 

eo = (1 + 6Ze)e, Aq = ^SZazZ + (1 + ^dZAA)A 
Th,o = Th + 6Th, Tho = Th + ST^ 



^H±,o = ^H±+'^^|r±, H+ = hzH+G+G+ + {1 + Uzh^h+)H+ (3.4) 

With the on-mass-shell conditions and tadpoles renormalization condition T = T + 5T = 0, we 
can obtain the renormalized constants expressed as, 

dZ^A = --Ke —2 |p2=o, o^^ZA = 2 2 ' (3.5) 



^Rv ^'W^R7 lz5 gSj^^(p') | , sin%i?e£|^ (0) 
2^^^^-^ + ^2^^^^ = 2^^^^'^^=° + cos^H- ml 



= --5Zaa + —T,^ZzA = -^Re — \p2=o + ^ ^2 > (3-6) 

5Zh+h+ = -Re^^\p2^M2^^, 5Mh± = ReT.f}\Ml^) - 5b aa, (3.7) 



where 



8b AA = 7, -[(JTff (sin^ /3 cos a-|-cos^ /3sina)-|-(5T/j(cos^ /3cos a— sin^ /3sin a)l, (3.. 

mw sm Ow sm 2p 



5Th = Th, 5Th = n. (3.9) 

The notation Re appearing in Eas. (|3.5|l . (|3.6|) and (|3.7jl . means taking the real part of the loop 
integrals appearing in the self-energy. 

We take the fine structure constant at the Z'^-pole as input parameter, Then we use the 



counter-term of the electric charge in DR scheme expressed as jl7| [TH| 

J . Ar/ „2 / A , i„„ <5M , AT.f „2 / A , i„„ <5 



6(47r)2 j '^'-^-^r \- ■ '"^ xj I ' ^^''^"^ ' m| 




(3.10) 

where we take Xf = niz when m/ < mz and xt = rat- e/ is the electric charge of (s)fermion and 
A = 2/e — 7 + log47r. N^j is color factor, which equal to 1 and 3 for (s)leptons and (s)quarks, 
respectively. 

The one-loop virtual corrections to 77 — > H^H^ is represented as 



1 pimax 

a^(J,77 - H+H-) = / di 2ReJ2m'')^M% (3.11) 



where tmax,min = i^H± ~ ^ I y ~ 4M^-|-S, and the summation with bar over head means 
the same operation as that appeared in Ea. (|2.5jl . is the renormalized amplitude for virtual 
one-loop corrections. After renormalization procedure, <7^ is UV-finite. Nevertheless, it still 
contains the soft IR singularities. The IR singularity in the is originated from the virtual 
photonic loop correction. It can be cancelled by the contribution of the real photon emission 
corrections. We shall discuss that in the following subsection. 

III. 2 Real Photon Emission Corrections 

We denote the real photon emission process as 

7(Pi) + 7(P2) ^ a^ivz) + ^~(P4) + 7(^), (3.12) 

where k = {k^,k) is the four-momentum of the radiated photon, and pi, p2, Ps and p^ are the 
four-momenta of two initial photons and final charged Higgs pair H^H~ , respectively. The real 




Figure 2: The real photon emission diagrams for the subprocess 77 H^H 7. 

photon emission Feynman diagrams for the process 77 — > H^H^^ are displayed in Fig|51 In our 
paper, we adopt the general phase-space-slicing method PU] to separate the soft photon emission 
singularity from the real photon emission process. By using this method, the bremsstrahlung 
phase space is divided into singular and non-singular regions. Then the correction of the real 
photon emission is broken down into corresponding soft and hard terms 

Affreai = ^^soft + ^^hard = ^oi^soft + ^hard)- (3.13) 

In the c.m.s. frame, the radiated photon energy = \J\k\^ + rni^ is called 'soft' if < AE^ 
or 'hard' if k^ > AE.y. Here, is a small photon mass, which is used to regulate the infrared 
divergences existing in the soft term. Although both A&soft and Aa^ard depend on the soft 
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photon cutoff A£'^/£'fe, where Ej, = ^ is the electron beam energy in the c.m.s. frame, the real 
correction A&reai is cutoff independent. In the calculation of soft term, we use the soft photon 
approximation. Since the diagrams in FigEJwith real photon radiation from the internal charge 
Higgs line or photon-charge Higgs vertex do not lead to IR-singularity, we can neglect them in the 
calculation of soft photon emission subprocesses (|3.12|) by using the soft photon approximation 
method. In this approach the contribution of the soft photon emission subprocess is expressed 

as mum 



dAasoft = -dao 



P3 Pi 



H 2 



(3.14) 



where the soft photon cutoff AE^ satisfies k'^ < AEry <C VS. The integral over the soft photon 
phase space has been implemented in Ref.|21j. then one can obtain the analytical result of the 
soft real photon emission correction to 77 H^H~ . 

As mentioned above, the IR divergence of the virtual photonic corrections can be exactly 
cancelled by that of soft real correction. Therefore, A(Tvii.-|_soft5 the sum of the virtual and soft 
contributions, is independent of the IR regulator m^. In the following numerical calculations, 
we have checked the cancellation of IR divergencies and verified that the total contributions of 
soft photon emission and the virtual corrections are numerically independent of m^. In addition, 
we present the numerical verification of that the total one-loop level EW correction to the cross 
section of 77 H^H^ , defined as Aa = Admr + AdreaU is independent of the cutoff AE^. 

Finally, we get an UV and IR finite correction Aa: 

Aa = Aa^r + Aa^-eai = f^oS 
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where 6 = 6mr + Sgoft + ^hard is the one-loop relative correction. 

III. 3 Calculation of the Parent Process e"'"e~ — >^ 77 ^ H^H~ 

The total cross section of the parent process e^e^ ^ 77 ^ H^H^ can be written as 

<t(s) = / dz— -^(t(77 ^ at s = z'^s) (3.15) 

ii?o/v^ 

with Eq = 2mu±, and being the e"'"e~(77) center-of-mass energy, ^^j^ is the distribu- 

tion function of photon luminosity, which is defined as: 



^§21 = 2^^"^ ^F,,,{x)F,i,{z'/x) (3.16) 

CLZ } -,'2 1 ^ X 



For the initial unpolarized electrons and laser photon beams, the energy spectrum of the back 
scattered photon is given by [23] 

^ 1 r 1 4x ^ ^ „^ 

^7/e = 7^[l-^+:^ Z-— - + —. —] (3-17) 



where 



L>(0 l-x C{l-x) C'^{l-x) 



where ^ = "le and -Eq are the incident electron mass and energy, respectively, ujq is 

the laser-photon energy, and x is the fraction of the energy of the incident electron carried 
by the backscattered photon. In our calculation, we choose loq such that it maximizes the 
backscattered photon energy without spoiling the luminosity via e~^e~ pair creation. Then we 
have ^ = 2(1 + V2), Xmax ^ 0.83, and D{i) = 1.8. 
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IV Numerical results and discussion 

We take the SM input parameters as nie = 0.511 MeV, ni^ = 105.66 MeV ^ 771^ = 1.777 GeV, 
mz = 91.188 GeV, mw = 80.425 GeV, m„ = 66 MeV, nic = 1.2 GeV, mt = 178.1 GeV, 
rud = 66 MeV, nis = 150 MeV, = 4.7 GeV, aew{m\)-^\-j^ = 127.918121]. There we use 
the effective values of the hght quark masses (m„ and md) which can reproduce the hadron 
contribution to the shift in the fine structure constant a^n„im\)\Zb\. 

The MSSM parameters are determined by using FormCalc package with following input 
parameters j26| : 

(1) The input parameters for the Higgs sector are the charged Higgs mass Mij± and tan/3. 
The output masses of other Higgs bosons are fixed by taking into account the significant ra- 
diative corrections (Actually the results are almost invariant quantitatively no matter which 
relation(tree or 2-loop level) we use, since those masse values of Higgs bosons are only adopted 
in the calculation of loop integrals.). 

(2) The input parameters for the chargino and neutralino sector are the gaugino mass param- 
eters Ml, M2 and the Higgsino-mass parameter /i. We adopt the grand unification theory(GUT) 
relation Mi = (5/3) tan^ 6'vi/-^2 for simplification |27j . 

(3) For the sfermion sector, we assume Mg = M^ = Mjj = M^ = M^ = MsusY and take 
the soft trilinear couplings for sfermions q and / being equal, i.e., Ag = Ai = Af. 

Except above SM and MSSM input parameters, we have to input some other parameters in 
our numerical calculations, for example, the IR regularization parameter and the soft cutoff 
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/^E^/E\,. In our following numerical calculations, we take AE^/Eb = 10~^ and rrij = 10 ^ GeV, 
if there is no other statement. As we know, the final results should be independent on the IR 
regulator m-y and the soft cutoff AEj/Ej,. For demonstration, we present the Fig.3, which 
shows the corrections to the cross section of the subprocess 77 — > H'^H~ versus the soft cutoff 
AEg/Ej, in conditions of M^± = 250 GeV, s/i = 1000 GeV and the input parameters in 
Set l(see below). The dashed, dotted and solid lines correspond to Aamr+soft, Aatard and the 
total one-loop electroweak correction Aa, respectively. As shown in this figure, the full one-loop 
EW correction Aa is independent of the soft cutoff AE^/E^ as AE-y/E^ running from 10~^ to 
10~^, although both Aa^ir^soft a^id Adhard depend on the soft cutoff strongly. 
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Figure 3: The full one-loop corrections to the subprocess 77 H^H~ as the functions of the 
soft cutoff AE/Eb. 

In order to show the numerical results and discuss the effects of the radiative corrections to 
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the subprocess of 77 — > H^H quantitatively, we choose the following four typical input data 
sets: 

Set 1: tan/? = 4, Mh± = 250 or 500 GeV, MsusY = 200 GeV, n = 600 GeV, M2 = 200 GeV 
and ^/ = 400 GeV. 

Set 2: tan P = 20, Mh± = 250 or 500 GeV, M^c/sy = 400 GeV, ^ = 1000 GeV, M2 = 200 GeV 
and ^/ = 800 GeV. 

Set 3: tan/? = 40, M^± = 250 or 500 GeV, MsusY = 200 GeV, /Li = 200 GeV, M2 = 1000 GeV 
and ^/ = 300 GeV. 

With the input parameters tan/3, Mfj±, MsusY-, A*, M2 and ^/ in one of the above data 
sets, we can obtain all the masses of super symmetric particles by using package FormCalcj26j. 
The input data 5'etl(or Set2) with small(or mediate) tan/?, makes the gaugino-like case with 
lighter(or heavier) sfermions, while the input data Set3 with larger tan/3 induces higgsino-like 
case. 

We also give the results at the SPSla' point from the SPA project j2Hl- The fundamental 
SUSY parameters in SPA project are compatible with all available precision data and actual 
mass and cosmological bounds. The SPA convention parameters are defined in the DR scheme 
at the scale of Q = 1 TeV. A translation from these parameters to our on- mass-shell definition 
can be performed by subtracting the corresponding counter terms, i.e. "pO^s _ "p^Q) — SV{Q). 
Then we get the pole mass of charged Higgs as Mfj± = 438.dGeV. For all other parameters 
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that do not enter in the tree-level calculations, either DR or OMS value can be used, since their 
difference is of higher order. 





Figure 4: The Born and the full one-loop level electroweak corrected cross sections of the 
subprocess 77 H^H~ versus c.m.s. energy \/I are plotted in Fig.4(a)(M^± = 250Gey) 
and Fig.4(c)(Mj:^± = bdQGeV)). The corresponding relative corrections as the functions of the 
c.m.s energy \/I are shown in Fig.4(b) and Fig.4(d), respectively. 

The Born and the full one-loop level electroweak corrected cross sections for the subprocess 
77 — > H^H~ as the functions of c.m.s. energy of 77 collision with above three input data sets are 



15 



displayed in Fig.4(a) with M^± = 250 GeV and in Fig.4(c) with Mjj± = 500 GeV, respectively. 
The corresponding relative corrections are depicted in Fig.4(b) and Fig.4(d). We can see that 
when s/l ~ 558 (1118) GeV, the tree level cross section reaches the maximal value 0.164 (0.041) 
pb in Fig.4(a) (Fig.4(c)). But its maximum value is shifted to 0.141 (0.337) pb after including 
the one-loop SUSY EW corrections. On the curve for the input data Set 1 in Fig.4(b) there 
exist small resonance spikes in the regions around the vicinities of VI ~ 2m(^ ~ 762.5 GeV and 
\/I ~ 2777,^2 ^ 968.3 GeV. On the curve for input data Set 2 the resonance peak is located at 
■\/§ ~ 2mj2 ~ 1137.1 GeV. For the curves of the input data Set 3, the resonance effect can be 
seen around the position of \/I ~ 2mf^ ~ 715.4 GeV. On the curves for the input data Set 2 
in both Fig.4(c) and Fig.4(d), we can see the resonance spikes at y/I ~ 2m|^2 ~ 1137.1 GeV. 
Both Fig. 4(b) and Fig. 4(d) show that the relative corrections have their maximal values at the 
position near the threshold energies and then decrease quantitatively with the increment of 
At the position of colliding energy V§ = 2 TeV shown in Fig.4(b) (Fig.4(d)), the relative 
electroweak correction can reach -34.5% (-34.8%), -39.4% (-56.4%) and -60.1% (-32.6%) 
for input data Set 1, Set 2 and Set 3, respectively. We can see from Fig.4(b) that the absolute 
relative corrections for the input data Set 3 can be rather large, and are generally larger than 
the corresponding ones for the input data Set 1 and Set 2. That is because in the input data 
Set 3 we have a very small sbottom mass m^^ = 76.67 GeV. While Fig.4(d) shows that the 
absolute relative corrections for the input data Set 2 are the largest among all the three input 
data sets, since there the conditions of mr fr f, ~ Mjj± = 500 GeVii = 1, 2) are satisfied. 

In Fig.5 we present the full one-loop relative electroweak corrections for the subprocess 
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Figure 5: The full one-loop relative electroweak corrections for the subprocess 77 — H'^H~ as 
the functions of the charged Higgs mass M^± . 

77 H^H~ as the functions of the charged Higgs mass M[j± with input data Set 1, Set 2 
and Set 3 separately. The corresponding collision energies in the c.m.s. of incoming 77, \/?, are 
1000 GeV, 1500 GeV and 2000 GeV ^ respectively. As shown in the figure all the curves are less 
sensitive to Mjj± , except in the region of Mjj± > 390 GeV. For the curve with the input data 
Set 3, the resonance effect can be seen around the positions of Mjj± ~ m^^ + m-^^ ~ 410.6 GeV 
and Mjj± ~ m^^ +"^61 ~ 426.2 GeV . In this figure the curve for the input data Set 1 shows that 
the relative correction is almost stable except in the energy region approaching the threshold 

~ 2Mfj± ~ 1000 GeV . However, the curve for the input data set Set 3 varies sharply in the 
region of 380 GeV < Mjj± < 440 GeV due to the resonance effects. 

In each figure of Fig. 6, Fig.7 and Fig.8, we take the input parameter sets as: Set 1, Set 2 and 
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Figure 6: The full one-loop relative electroweak corrections for the subprocess 77 H'^H as 
the functions of tan/3. 

Set 3, except the charge Higgs mass being 200 GeV, 300 GeV and 500 GeV, the colHding energy 
being 500 GeV, 1000 GeV and 2000 GeV, respectively. In Fig.6 we present the full one-loop 
relative electroweak corrections for the subprocess 77 H^H~ as the functions of the ratio of 
the vacuum expectation values tan /3. We can see from the curves for input data Set 1 and Set 2 
that the relative corrections decrease slowly with the increment of tan/3 except the curve for 
Set 2 in the region of tan j3 > 32. For the curve of the input data Set 3, the resonance effect can 
be seen around the position of tan /3 ~ 21 where the condition of Mfj± ~ m^^ + m^^^ ~ 500 GeV 
is satisfied. 

In Fig. 7 we present the full one- loop relative electroweak corrections for the subprocess 
77 — H^H~ as the functions of MguSY- We can see that the relative corrections are not 
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Figure 7: The full one-loop relative electroweak corrections for the subprocess 77 H'^H~ as 
the functions of MsuSY- 

sensitive to MsuSY except at the vicinities of the resonance points. For the curve of the input 
data Set 1, the resonance effect can be seen around the position of MguSY ~ 276 GeV where 
a/I ~ 2mj^ ~ 500 GeV . For the curve of the input data Set 2, the resonance effect can 
be seen around the positions of MguSY ~ 346 GeV (where Mjj± ~ mj^ + m^^ ~ 300 GeV) 
and MsusY ~ 594 GeV (where \f§ ~ 2m(^ ~ 1000 GeV). For the curve of the input data 
Set 3, the resonance effect can be seen around the positions of MguSY ~ 230 GeV (where 
Mh± ~ + ~ 500 GeV) and MsuSY ~ 242 GeV (where M^± ~ m^-^ + m^;^ ~ 500 GeV). 

In Fig.8 we present the full one-loop relative electroweak corrections to the subprocess 77 — 
H'^H" as the functions of M2. We can see that the relative corrections are not sensitive to 
M2 except the curve for the input data Set 3. On the curve of the input data Set 3, there are 
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Figure 8: The full one-loop relative electroweak corrections for the subprocess 77 H'^H~ as 
the functions of M2. 

five resonance points which come from the fact that the masses of some chargino and neutralino 
are lighter than Mh±. The resonance effect can be seen around the positions of M2 ~ 244 
GeV (where Mff± ~ m-+ + m^o ~ 500 GeV), M2 ~ 280 GeV (where Mh± ~ m-+ + m^o 
500 GeV), M2 ~ 288 GeV (where Mh± ~ m~+ + m^o ~ 500 GeV), M2 ~ 334 GeV (where 
Mh± ~ m-+ + m^o ~ 500 GeV) and M2 ~ 620 GeV (where Mh± ~ m~+ + m^o ~ 500 GeV). 

A,2 ''-1 A,l '^3 

For the input data set Set 1 and Set 2, the relative corrections are almost stable, the relative 
corrections are about —16.9% and —27.6%, respectively. 

Fig.9(a) is the plot of the Born and the full one-loop level electroweak corrected cross sections 
for the parent process e'^e" — >^ 77 — H^H~ versus the electron-positron colliding energy with 
Mh± = 250Gey. At the position of y/s ~ 1.5 TeV the cross sections reach their maximal values. 



set1 
set2 
sets 
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Figure 9: The Born and the full one-loop level electroweak corrected cross sections for the parent 
process e+e" -^■ 77 — > if+if^, are shown in Fig.9(a)(Mj:^± = 25QGeV) and Fig.9(c)(Mj^± = 
500Gey), respectively. The corresponding relative corrections as the functions of the c.m.s 
energy of the incoming electron-positron pair, are shown in Fig.9(b) and Fig.9(d), separately. 
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e.g., for the curve of the input data Set 1, the cross section has the maximal value 52.2 fb. When 
y/s > 1.5 TeV, the cross sections decrease slowly with the increment of the colliding e+e" c.m.s 
energy. Fig.9(b) shows the corresponding relative corrections as the functions of colliding e+e~ 
energy. We can see that the absolute relative corrections increase obviously with the increment 
of s/s. When \/i = 2 TeV, the relative correction can reach its maximal value —29.2% for 
the curve of the input data Set 3. Fig.9(c) is the plot of the Born and the full one-loop level 
electroweak corrected cross sections for the parent process e+e~ — 77 — s- H^H~ versus the 
electron-positron colliding energy with M^± = SOOGeF. The cross sections increase with the 
increment of the colliding c.m.s energy. When ^/s = 2 TeV , the corrected cross section reaches 
its maximal value 10.9 fh for the input data Set 3. Fig.9(d) shows the corresponding relative 
corrections as the functions of the colliding e"^e~ energy. We can see that the absolute relative 
corrections increase apparently with the increment of -y/s. When -y/s = 2 TeV, the relative 
correction can reach —47.3% for the input data Set 3. 

Fig. 10 (a) and Fig. 10(b) shows the Born and the full one- loop level electroweak corrected 
cross sections and the corresponding relative corrections for the subprocess 77 H^H~ as the 
functions of the c.m.s energy a/I at the SPSla' point. From Fig.lO(a) we can see that the cross 
sections decrease with the increment of the colliding c.m.s energy \f§ when \/i > 1 TeV. At 
the point of \/I = 980 GeV , the tree level and one-loop level corrected cross sections are 53.2 fh 
and 49.4 fh, respectively. When ^/I = 2 TeV, the tree level and one-loop level corrected cross 
sections go down to 21.0 fh and 15.7 fh, respectively. In Fig.lO(b) we can see that there are 
two small peaks at the vicinities of ~ 2mc2 ~ 1048.3 GeV and \/i ~ 2m( ~ 1137.4 GeV . 
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Figure 10: The Born and the full one-loop level electroweak corrected cross sections for the 
subprocess 77 — > i^+i?" (parent process e+e" — > 77 ^ H~^H~), and their corresponding 
relative corrections as the functions of the c.m.s colliding energy V§{y/s) at the SPSla' point 
are shown in Fig.lO(a) and Fig. 10(b) (Fig. 10(c) and Fig.lO(d)), respectively. 
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The absolute relative correction increases with the increment of v^. When goes from 1 TeV 
to 2 TeV, the relative correction varies from —7.76% to —19.5%. Fig. 10(c) and Fig. 10(d) are 
the plots for the cross sections and the corresponding relative corrections of the parent process 
e'^e~ — > 77 ^ H~^H~ as the functions of the c.m.s energy of the incoming electron-positron 
pair, separately. Prom Fig.lO(c) we can see that the cross sections increase with the increment 
of the colliding e+e~ c.m.s energy y/s. When a/s = 2 TeV, the tree level and one-loop level 
corrected cross sections are 18.7 fb and 16.8 fb, respectively. In Fig.lO(d) we can see that the 
absolute relative correction increases with the increment of i/s. When y/s goes from 1.1 TeV 
to 2 TeV, the relative correction varies from —0.86% to —10.2%. 

V Summary 

In this paper, we present the calculation of the full one-loop electroweak corrections to the 
subprocess 77 H^H~ and parent process e+e~ — >^ 77 — >^ H~^H~ at a linear collider in the 
MSSM. We analyze the dependence of the relative corrections for the subprocess on colliding 
energy, charged Higgs boson mass and several supersymmetric parameters. We find that these 
corrections generally reduce the Born cross sections and the relative corrections are typically 
few dozen percent for both the subprocess and the parent process . With the input data Set 3, 
the relative corrections to the subprocess are obviously sensitive to Mh± , tan (5, MsuSY and M2 
in some parameter space due to the resonance effects. However, with the input data Set 1 and 
Set 2, the relative corrections to the subprocess are less sensitive to these parameters comparing 
with the curves with input data Set 3. We also give the numerical results at the SPSla' point, 
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it shows that with ^ varying from 1 TeV to 2 TeV, the relative correction to the subprocess 
runs from —7.76% to —19.5%. We conclude that the complete one-loop electroweak corrections 
to both subprocess 77 — H^H~ and the parent process e^e~ — 77 — H^H~ are generally 
significant and should be considered in the precise analysis. 
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Figure Captions 

Figure 1 The leading order diagrams for the 77 — > H'^H~ subprocess. 
Figure 2 The real photon emission diagrams for the subprocess 77 H^H^'j. 
Figure 3 The full one-loop corrections to the subprocess 77 — > H^H^ as the functions of 
the soft cutoff A£;/^fc. 

Figure 4 The Born and the full one-loop level electroweak corrected cross sections of the 
subprocess 77 H^H^ versus c.m.s. energy are plotted in Fig. 4(a) (Mj:^± = 250GeV) and 
Fig. 4(c) (Mj:^± = 500GeV)). The corresponding relative corrections as the functions of the c.m.s 
energy \/I are shown in Fig. 4(b) and Fig. 4(d), respectively. 

Figure 5 The full one-loop relative electroweak corrections for the subprocess 77 H^H~ 
as the functions of the charged Higgs mass Mjj± . 

Figure 6 The full one-loop relative electroweak corrections for the subprocess 77 H^H~ 
as the functions of tan /3. 
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Figure 7 The full one-loop relative electroweak corrections for the subprocess 77 H^H~ 
as the functions of MsuSY- 

Figure 8 The full one-loop relative electroweak corrections for the subprocess 77 H^H~ 
as the functions of M2. 

Figure 9 The Born and the full one-loop level electroweak corrected cross sections for 
the parent process e^e~ — ^ 77 — > H'^H~, are shown in Fig.9(a)(Mj:^± = 250GeV) and 
Fig.9(c)(M^± = bOOGeV), respectively. The corresponding relative corrections as the func- 
tions of the c.m.s energy of the incoming electron-positron pair, are shown in Fig.9(b) and 
Fig.9(d), separately. 

Figure 10 The Born and the full one-loop level electroweak corrected cross sections for 

the subprocess 77 — > H~ {parent process e+e~ — >• 77 ^ H~^H~), and their corresponding 
relative corrections as the functions of the c.m.s colliding energy ■\/§{^/s) at the SPSla' point 
are shown in Fig.lO(a) and Fig.lO(b)(Fig.lO(c) and Fig. 10(d)), respectively. 
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